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Abstract 

We develop a novel approach to chiral meson-baryon dynamics incorporat- 
ing unitarity constraints and explicit resonance fields. It is based on the most 
general structure of any pion-nucleon partial wave amplitude neglecting the un- 
physical cuts as derived from the N/D method. This amplitude is then matched 
to the one-loop heavy baryon chiral perturbation theory result at third order 
and to tree level exchanges of baryon- and meson states in the s, t and u chan- 
nels. This generates the left-hand cuts. The unitarization procedure does not 
involve form factors or regulator functions. The resonance parameters are de- 
termined from fits to the S- and P-wave pion-nucleon partial wave amplitudes 
for energies up to 1.3 GeV. In particular, the A(1232) is accurately reproduced 
whereas scalar and vector meson couplings are less precisely pinned down. We 
also obtain a satisfactory description of the 7rN threshold parameters. Further 
extensions of this method to coupled channels and the three-flavor case are 
briefly discussed. 
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1 Introduction 



There are many reasons why it is interesting to study pion-nucleon scattering. First, in the 
threshold region, chiral perturbation theory is applicable and thus the chiral structure of 
QCD can be investigated, see e.g. [jjj (and references therein). This very precise method 
is, however, only applicable in the threshold and in the unphysical region since there the 
pion momenta are small. At higher energies, additional physics due to the appearance of 
resonances and coupled channel effects sets in. In fact, the dispersive analysis of pion- 
nucleon scattering data has been and still is the best method to analyze the properties of 
the baryon resonances (as long as they couple strongly to the 7rN system). Many models 
and approaches of different types have been developed to deal with these phenomena, we 
just mention the meson-exchange || and the extended tree level [||] models. While all these 
models in general give a good description of the phase shifts (and inelasticities), they are 
not including the pion loops resulting from the chiral structure of QCD in a systematic 
fashion (the often employed unitarization schemes or resummation techniques include some 
classes of these diagrams, but not all). Furthermore, in these references most of the tree 
level diagrams involving the exchange of meson resonances do not have the momentum 
dependence required by chiral symmetry. Consequently, it would be interesting to devise a 
scheme that at low energies exactly reproduces the chiral perturbation theory amplitudes 
(to a given order in the chiral expansion) but also allows one to work in the resonance region, 
including also resonance exchange at tree level in conformity with the non-linearly realized 
chiral symmetry. Such an approach will be developed here. It is similar to the one set up 
in ref.|| for meson-meson scattering and will be explained in more detail below. Before 
doing that, it is worth to stress that this method can easily be extended to the coupled 
channel case and also to the three-flavor sector, in particular also for the description of 
kaon-nucleon scattering. While there have been already some successful studies of chiral 
SU(3) dynamics with coupled channels for this and related processes || 0], no attempt has 
been made so far to match the 7rN subprocesses to the accurate chiral perturbation theory 
predictions. Thus, the model presented here can be considered as a first building block for 
setting up a relativistic coupled channel scheme which fulfills these requirements. 

It is known that dispersion theory can be used to match on chiral perturbation theory and 
thus allows one to extend the range of energies one can consider. Of course, any unitarization 
scheme is at variance with the power counting underlying an effective field theory like e.g. 
CHPT. It is, however, not known how to extend such a concept to include resonances (with 
the exception of the A(1232) isobar, see ref.||). On the other hand, unitarization methods 
allow one to exactly fulfill unitarity (which is only respected perturbatively in CHPT) and 
can generate S-matrix poles corresponding to resonances. We attempt to combine these 
concepts in a way that makes best use of these various features. For that, employing the 
well-known N/D method |J, we derive the most general structure for any 7rN partial wave 
amplitude when the unphysical cuts are neglected. Then, similarly as in ref. the 
unphysical cuts will be treated as a perturbation in a loop expansion with the HBCHPT 
power counting. Hence, our starting point of no unphysical cuts at all is the zeroth order 
approach in this scheme from which the first order approach, the one necessary to match 
with 0(p 3 ) HBCHPT, is developed. Note that tree level dia grams involving local operators 
and particle exchanges in all the channels are also included since they do not involve loops. 
This should be a good starting point for the channels which are dominated by unitarity and 
the presence of s-, u— or i-channel resonances, like e.g. the P33 partial wave which is largely 
given by the A(1232). The argument is further strengthened by the observation made in 
ref. (ll]], namely that the dimension two low-energy constants of the chiral effective 7rN 
Lagrangian are saturated by s— , u— and t-channel resonance contributions, most notably 
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the A and the p(770). In addition, for the scalar-isoscalar channel one can deduce a 
certain contribution from a meson with these quantum numbers. In ref.[ll], this meson 
was taken to be a light a, but in fact only the ratio of the scalar meson-nucleon coupling 
constant to the scalar meson mass plays a role. We will come back to this point in much 
more detail below. In addition, we will also include the Roper iV*(1440) resonances which 
plays a minor but non-negligible role in the energy range we will consider (since we do 
not calculate inelasticities, we stay below center-of-mass energies of about 1.3 GeV). Note 
that preferably we should match to a relativistic 7rN chiral perturbation theory amplitude 
because in that way one can avoid the difficulties with the HBCHPT approach as detailed 
in ref . ] 12 ] - Such an amplitude, based on a proper power counting, is not yet available but 
under construction, see ref. [|l2j . It remains to be seen whether the matching to such an 
amplitude will lead to better results than the ones presented here. 

Our manuscript is organized as follows. In section 2, we collect the formalism necessary 
to describe elastic pion-nucleon scattering. Section 3 we construct the basic tree level 
diagrams supplemented by the 0(p 3 ) HBCHPT loops. In section 4 the unitarization method 
is discussed and in sections 5 the results are presented. The most relevant conclusions of 
the work are summarized in section 6. 



2 Prelude: Kinematics 

In this section, we collect the necessary formalism concerning elastic pion-nucleon scat- 
tering. This serves to set our notation and to keep the manuscript self-contained. In the 
following, we consider the 7rN scattering amplitude. In the center-of-mass frame (cm.), the 
amplitude for the process n(qi) + N(q2) — > tt((/2) + Nfa) can be written as (suppressing 
isospin indices) 



T(W,t) 



E + m 
2m 



g(W,t)+ia-(q 2 x q\)h(W,t) 



(2.1) 



where q\ and q\ are the three-momenta of the incoming and the outgoing pion, respectively, 
W = y/s is the total cm. energy, m is the nucleon mass and t = (q± — q2) 2 is the invariant 
momentum transfer squared. Furthermore, g(W,t) refers to the non-spin-flip amplitude and 
h(W, t) to the spin-flip one. For a given total isospin (J) of the 7rN system, the partial wave 
amplitudes T^, where I refers to the orbital angular momentum and the subscript ± to the 
total angular momentum ( J = I db 1), are given in terms of the functions g(W, t) and h(W, t) 
via 

TL(W) = ^±^ J^dzlgP^ + lqfhiP^-zPtiz))] , (2.2) 

where z = cos#, with 9 the cm. angle between the incoming and the outgoing pion, \q\ 
is the modulus of the cm. three-momentum of any of the pions, and E = \J m 2 + q 2 is 
the cm. nucleon energy. Note that s, t and u are the Mandelstam variables subject to the 
constraint s + t + u = 2m 2 + 2M 2 , with M the pion mass. 

As a consequence of unitarity, any partial wave amplitude can also be written in terms of 
a corresponding phase shift 6f ± as 

TL(W) = ^ exp(Wk) sin 5/ ± . (2.3) 
This implies the unitarity relation 

ImT/ ± (H/) = ^L|T/ ± (^)| 2 , (2.4) 
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which is valid above threshold and before the opening of any inelastic channel. 

The amplitudes g(W, t) and h(W, t) are specially suited for a non-relativistic treatment of 
the 7rN scattering, as it is carried out e.g. in Heavy Baryon Chiral Perturbation Theory 
(HBCHPT) (for a review on HBCHPT see [13| and an update is given in [|l4]]).[^] However, 



in a fully relativistic approach of 7rN scattering, it is more convenient to work with the 
invariant amplitudes A(s,t) and B(s,t), which are free of kinematical singularities and are 
more suited for studying the analytical properties of the scattering amplitude Q . In terms 
of the invariant functions A(s,t) and B(s,t), the 7rN scattering amplitudes are written as: 



T(W,t) =u(;P2,A 2 ) 



A(s,t) + -(q 2 + qx)^B( S ,t) 



u(pi,Ai), (2.5) 



where the Aj, appearing in the Dirac spinors, denote the helicities of the incoming/outgoing 
nucleon. One can derive in a straightforward way the relation between the A, B and g, h 



amplitudes. From the expressions given in ref. [15] one obtains 



g = — — - — -75- U4m(E + m) -t) A+ \(W + m)(t + 4|<f| 2 ) + 4mw(E + m)] B\ , 
2(m + Ey L L 

h = ( ]^ W + m ) B - A }^ ( 2 - 6 ) 
(m + Ey 

with w the cm. pion energy, w = W — E. It is evident from the previous equations 
that, since the amplitudes A(s,t) and B(s,t) are free of kinematical singularities, the same 
does not hold for the amplitudes g(W, t) and h(W, t) due to the dependence on W = 
y/s. Note that the cut singularities induced by this dependence have nothing to do with 
the exchange of any particle states. This is the reason why they are called kinematical 
singularities, in contrast to those due to the exchange of physical particles, which are called 
dynamical singularities. In the invariant amplitudes A(s, t) and B(s, t), only this latter type 
of singularity appears. 

So far, we have ignored the isospin of the pion (I n = 1) and the nucleon (Ijy = 1/2) fields. 
Consequently, in 7rN scattering there are two isospin amplitudes, for total isospin 1=3/2 
and 1=1/2, respectively. These amplitudes can be written in terms of the ir + p — * tr + p 
amplitude, Ti, and of the one for n~p — ► ir~p, ?2, as: 

T I=3/2 (W,t) = T x (W,t) , (2.7) 
T I=1/2 (W,t) = ^T 2 (W,t)-~T 1 (W,t) . 

If we denote by Ai(s,t) and Bi(s,t) the invariant amplitudes for T±(W, t) and by A2(s,t) 
and B2(s,t) the corresponding ones for T2(s,t) then, by crossing, the following relations 
arise: 

A 2 (s,t) = +A 1 (u,t) , (2.8) 
B 2 {s,t) = -Bi(u,t); 

which allow one to consider just the ir + p — > ir + p reaction in order to determine the 1 = 1/2 
and 3/2 scattering amplitudes. This is similar to the case of elastic pion-kaon scattering, 



see e.g. ref. 1 16] 



# 3 We remind the reader on the remarks made in the introduction on the calculation of pion-nucleon 
scattering in a relativistic formulation of baryon chiral perturbation theory. 
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3 Tree level amplitudes and chiral loops 



In this section we construct the basic amplitudes which we will unitarize in the next section. 
First, we consider the tree level contributions to T\ arising from the lowest order meson- 
baryon Chiral Perturbation Theory (CHPT) Lagrangian, jC^jv - Then, the exchange of the 
delta isobar A(1232) and of the iV*(1440) and of meson resonances in the t-channel is 
also taken into account. Note that beyond leading order, the chiral Lagrangian contains 
local contact terms of second (and higher) order. In our approach, these will be generated 
by the explicit contribution from the aforementioned resonances. That this is a sensible 
starting point is based on the observations made in ref.[ll], where the dimension two low- 
energy constants were fixed from data and it was shown that resonance saturation allows 
to explain these values. Having set up the tree level Lagrangian, we will then consider the 
0(p 3 ) HBCHPT loop contributions [L7|, Q to the different isospin amplitudes taking the 
results directly from ref.O. 



3.1 Pion— nucleon Born terms 

In the two-flavor SU(2) case, the one we are interested in here, the lowest order CHPT 
pion-nucleon Lagrangian takes the form: 

4iv = ~m+ ^ 7 M 75^)^ (3-1) 

In this expression one has to take into account that the nucleon mass, m, and the axial- 
vector coupling gA are to be taken in the chiral limit (m u = = 0, m s fixed). However, 
in the following we will use the physical values of these quantities since, when including the 
chiral loops at 0(p 3 ) j!]], they will be renormalized, up to this order, to their physical values. 
In the following we will take gA = 1-26 and m = 938.27 MeV. The choice of the proton mass 
is justified because we are only considering the scattering of pions off the proton. From this 
process we derive the two isospin amplitudes 7rN scattering, as explained above. On the 
other hand, the proton p and the neutron n fields are collected in the isospinor ^ 

•-(;), 

the covariant derivative acting on the nucleon field is given by (note that we do not 
consider external fields and isospin breaking here) 

= + , 

r M = \[u\d^u] . (3.3) 

Here, Un = i u^dpUu^ and U = u 2 = exp[fr ? 7T ? /.F 7r ], where t j refers to the Pauli matrices (in 
isospin space). The tree level diagrams from C^ N are depicted in fig.l. The corresponding 
expressions are: 

&r£ \ u — m z J 

t\ = k 



4FJ 



with F n = 92.4 MeV the (weak) decay constant of the pion and i? M = (q± + q^)^- Note that 
the direct nucleon exchange diagram does not contribute to 7r + p — > n + p. 
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n 



a) b) 

Figure 1: Diagrams from C^ N : l.a) Crossed nucleoli exchange, l.b) Seagull. 
3.2 Baryon resonances 

Next, we consider resonance excitations in the s- (or it-)channel. We treat explicitly the 
spin-3/2 A(1232) and also the first even-parity excited state of the nucleon, the Roper 
iV*(1440). Most important for the energy-range we are considering is the A(1232) reso- 
nance. It contributes to 7rN scattering through the A(1232) exchange diagrams, see fig. 2. 
The interactions of the A field with pions and nucleons are given by the following effective 
Lagrangian, which is the leading term of an appropriate chiral invariant interaction [13] 

CattN = ^A^e uX (Z)Nd x n + h.c. , (3.5) 



where g& n N is the AvriV coupling and M = 139.57 is the charged pion mass. This choice is 
justified since we are considering elastic ir + p scattering. On the other hand, T is the \ — ► | 
isospin transition operator which satisfies 

3 3 1 

l^T^-X^^-X^Tl = 5ab ~ 3 TbTa ' ^'^ 
Furthermore, the Dirac matrix operator 0„„(Z) is given as 



9^(Z) = g„ u + Q(l + AZ)A + z\ lixlv = - (Z + , 



(3.7) 



with the right-hand side valid for A = — 1. The quantity Z is called an off-shell param- 
eter, which enters the fully relativistic AttN vertices. Finally, the spin-3/2 propagator 
(momentum P, v to \i) is: 

" ~ 3 7 ^ ~ + 3m A ) ' (3 ' 8) 

with m A the mass of the delta resonance, we take m A = 1232 MeV. Note that this mass 
used in the tree graphs is not necessarily equal to the real part of the S-matrix pole. We 
will come back to this point later on. For the diagram represented in fig. 2a, corresponding 
to the direct exchange of the A, we find 



Tf- a = Af a + B^ a , 

A t M = r — ^wt~ iTT W 2 ™ [6M 2 - 2s - 3t - Ask 2 + m 2 (2 + 4k 2 )] 

0771a M z [s — m A ) 

+ m A 3 [6M 2 - 2s - 3t + Ask - 8sk 2 + 2m 2 (1-2k + 4k 2 )] 

+ 2 m A [-M 4 - M 2 s - 2 s 2 k + 4 s 2 k 2 + m 2 (M 2 + 2 s k - 4 s k 2 )] 

- m [m 4 + M 4 + 2 M 2 s + s 2 - 4 s 2 k 2 - 2 m 2 (M 2 + s - 2 s k 2 )] } , 
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a) b) 

Figure 2: Diagrams with the exchange of the A(1232) resonance. 

B A.a = 9A^ f m 4_( M 2 + s _ 2sK ) 2 + 2m ?(M 2 + S-2sK-2 S K 2 ) 

6mA M z (s — m A ) L 7 

+ 4 m A 3 m (l - 2 k + 4 k 2 ) + 2 m A m (m 2 -M 2 -s + 4sk-8sk 2 ) 

+ m A 2 [4M 2 -3t-4M 2 k-4sk + 4sk 2 + 4m 2 (1 + k + k 2 )]} , (3.9) 

with k = Z + ^. The diagram of fig. 2b can we obtained from the former one just by crossing. 
The result is: 

A?- b (s,t) = +±A?- a (u,t), 

Bf- b (s,t) = -±B?- a (u,t). (3.10) 

Let us consider now the exchange of the iV*(1440) Roper resonance, cf. fig. 3. Although 
this resonance is heavier than the A(1232) and not as strongly coupled to the 7rN system, 
in ref.[^] it was shown that the Roper octet can have an appreciable interference with 
the contributions coming from the exchange of scalar resonances. The Lagrangian for the 
coupling iV*iV7r is Q 

C N *N* = ^-VRVn^-Ysu^n + h.c. (3.11) 



with VR = 0.53 ± 0.04 j2l| from the width of the TV* (1440) determined in ref.[g2| (note 
that the often used Breit-Wigner fits to this resonance lead to an appreciably larger width) . 
The resulting amplitude is given by: 

T l* = w2i AR 2^ \{u-m 2 ){m + m R ) + ^{u + m 2 + 2mm R )\ (3.12) 
8F4(u-m R ) { 2 J 

where m R = 1440 MeV is the mass of the Roper resonance. 



71 - 7 71 



p N* P 

Figure 3: s-channel exchange of the Roper resonance. 
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3.3 Meson resonances 



We now turn to the i-channel exchange of meson resonances. First, we consider the ex- 
change of an isosinglet scalar resonance, cf. fig. 4. This state should account for both the 
contribution of the lightest scalar singlet and octet 1 = states, which in large N c (where 
N c denotes the number of colors) should be degenerate. The SU(2) Sirir interaction can be 
written as flTlll 



£-Sn7T = S [c rn Tr(x+) + c d Tr^i^)] 



(3.13) 



with x+ = diag(M 2 , M 2 ) and the Cd, c m are two lowest order coupling constants. In ref.||, 
where a simultaneous fit to the all SU(3) related S-wave meson-meson scalar channels was 
done combining CHPT with the exchanges of resonances in the s-channel via the N/D 
method, the scalar singlet appears with a mass around 1 GeV and the octet is higher in 



energy, around 1.35 GeV. In ref. [25| it was shown that the requirement that the tree level 
Kit, Kr) and K rf scalar form factors vanish at infinity fixes these coupling in terms of the 
pion decay constant, 

p 

c~d = Cm = (3.14) 
V6 

Requiring also the saturation of the 0(p 4 ) CHPT counterterms for the former couplings, 
the nonet mass is then fixed to be around 1.2 GeV, which lies in an intermediate value 
between the singlet and octet mass given in ref.Pj. The well-known cr(550) meson, see e.g. 
refs. |^6L [5|, B7| should be considered as a reflection of the strong S-wave -kit interactions, loop 
physics and hence subleading in large N c counting rules (this has already been stressed a long 
time ago, see p8| ). Other QCD inspired approaches, which also establish that the lightest 
preexisting scalar meson should be around 1 GeV, are collected in refs.p9|) 3Cj. While the 
contributions of a preexisting scalar resonance appear at 0(p 2 ) in the chiral counting, the 
crossed pions loops, and hence their associated a pole, begins at 0(p 3 ). The coupling of 
the scalar 1 = resonance to the nucleon is calculated easily from the Lagrangian 

Csnn = -gsV* ■ (3.15) 
From eqs.( p.l3 ),( 3.15| ) we obtain the scalar meson exchange contribution to 7rN scattering, 

c m 2M 2 +c d (t-2M 2 ) 



T? = -2g s 



F 2 (t-M 2 



(3.16) 



with Ms the mass of the scalar resonance. 
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P P 

Figure 4: t-channel exchange of an isosinglet resonance, S. 



Next, we consider the exchange of the p(770) meson in the t-channel, fig. 5. Following 
refs. [18, 19], we consider the tensor representation for the vector fields. In this way, we 



generate terms in the chiral expansion of the pion-nucleon scattering amplitude that begin 
at 0(p 2 ). On the other hand, if we were to use the vector representation, the contribu- 
tions from the exchange of the p would begin to appear only for orders higher than 0(p 2 ). 
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Consequently, extra local contact terms would be needed in order to make both representa- 
tions, vector and tensor, equivalent [23, 19]. The tensor field representation is most useful 
in constructing chiral invariant building blocks for spin-1 mesons coupled to nucleons, pi- 
ons and external fields (like e.g. photons). The pir +/ n~ vertex can be obtained from the 



Lagrangian [18] 



(3.17) 



where Wau refers to the tensor field representing the octet of the lowest-lying vector reso- 
nances, 

/ „° . ,„„ 

P + 



( P° I ^8 

V2 ^ V6 



w, 



[IV 



p 

K* 



£ i_ o>8 



K*+ \ 
K*° 

2^8 
V6 / 



such that 



(0\W^\V,p) = iMy l [p^ u (p) -p,e»] 



(3.18) 



(3.19) 



with ea(p) the usual polarization vector for the vector state \V,p) with a mass in the chiral 
limit of My ~ M p = 770 MeV [pi] . The propagator is given by 



Ml 



Ml 



p z 



l€ 



[9n P 9va(M v - p 2 ) + g^PuPa - g^PvPp - v )] ■ 



(3.20) 



The coupling Gy can be obtained from the decay p — > 2tt and one obtains Gy = 69 MeV. 
Another possible way of determining Gy is from fitting the charge radius of the pion. That 
leads to a reduced value of Gy = 53 MeV [24]. We will distinguish below between both 
values when we present our results. One has to take into account that in ref. []l8|, |l^] 
an SU(3) notation is considered. We will also use this three-flavor notation to work out 
the p-couplings. The only difference for our considerations is that the Pauli matrices, r J , 
are changed to the Gell-Mann matrices A*. In this way the SU(3) representation of the 
Goldstone boson octet reads 



U 



exp i 



/ 7T° I VB 



7T 



IT 

K~ 



+ 



i la 



V2 Vq 
K° 



K+ 

V6 



(3.21) 



7t + , 
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P P 

Figure 5: t-channel exchange of the p(770) resonance. 



On the other hand, the pNN vertex can be inferred from ref. [IS] from the general vector 
to baryon octet coupling via 

Cvnn = R D / F Tr[Ba^(W^B))+S D/F Tr[B^(d u Wp P ,B)) 

+ T D/F Tr[B^(d x W^d x d v B)) + U D/F Tr[Ba x "(W^,d x d^B)] . (3.22) 
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Here, the symbol (A, B) denotes either the commutator, (A,B) = [A, B], or the anticom- 
mutator, (A, B) = {A, B}. In the first case, the coupling constant has the subscript F, in 
the second one the subscript D applies. For instance, the first term in the previous equation 
thus means 

R D Tr[Ba^{W^,B}]+R F Tr[Ba^[W pv ,B]} . (3.23) 

From the Lagrangians given in eqs.( |3.17D and ( 3.22j ), one obtains for the diagram depicted 
in fig. 5: 



2G V V2 



Ra 



2mR + + 



m U-l 



t 2 T+ 



(3.24) 



with R+ = Rp + Rp and analogously for S+, U+ and T + . All the former unknown param- 
eters, except T + , can be recast in terms of more familiar ones, by comparing eq.( |3.24 ) with 
the result that one would obtain by making use of the standard vector representation for 
the p meson. In this case, instead of the Lagrangian given in eq.(|3.22|), one has: 



C-pNN = ~ QpNN * 



2m . 



(3.25) 



with g p NN the vector coupling of the p to two nucleons and n p the ratio between the tensor 
and the vector couplings. After these preliminaries, the final expression that we obtain for 
the p-exchange contribution is given by: 



t Gy g p NN 

M p F%(t-Mj) 

QpNN K p Gy 



FlM n 



Vi u — s 

T7 + -a 

2 Am 



+ 



4m 

t 2 V2G v 
AF*{t-Mj) 



(3.26) 



Note that the coupling ~ t 2 T + has not been considered before although it is solely based 
on chiral symmetry. Its relevance for the pNN phenomenology will be discussed below. 



3.4 Chiral pion loops 

As explained in detail in the introduction, we consider here the 0(p 3 ) one-loop contributions 
obtained in the framework of heavy baryon chiral perturbation theory. After renormalizing 
all the bare parameters, these loop contributions are given in ref. Their contributions to 
the spin-non-flip and spin-flip amplitudes with given total isospin (I = 3/2 or 1/2), 53/2, 
g 1/2 and h 3/2 , h 1/2 are: 

53/2 loop = 
#1/2 loop = 

where we have followed the notation of ref. (l| to denote their calculated loop functions. 
For the h 3/2 ioo P and h i/2 loop functions analogous relations with the functions ^ op , ^i oop 
as the ones shown in eq.( |3.27j ) hold. 



^loop 



9i, 



oop 



^loop 

+ 2^ 



oop 



(3.27) 



4 Unitarization 



As a consequence of eqs^^) and ([2.6|), any 7rN partial wave amplitude has besides the 
dynamical cuts, which are due to the exchange of particle states in the s— , t- or u— channel, 
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further kinematical singularities arising from the y^-dependence in eq. (|2.6| ). In the follow- 
ing, we will consider a partial wave amplitude as a function of W, instead of s, and hence 
the kinematical singularities disappear. The dynamical cuts are usually divided into two 
classes: the physical or unitarity cut and the unphysical cuts. The physical cut is required 
by unitarity eq. fl2.4| ), and the unphysical cuts are due to the exchanges of particles in the 
crossed channels. A more detailed discussion of the ideas underlying this approach are given 
in ref.Q. Here, we are concerned with the extension of such a scheme to the pion-nucleon 
sector. 



4.1 The right— hand cut 

We consider first the physical or right-hand cut and neglect, for the moment, the unphysical 
ones. In order to obtain a unitarized 7rN amplitude, we follow ref.|J, in which meson-meson 
scattering was considered. In this work, making use of the N/D method ||, the most 
general structure of a meson-meson partial wave amplitude, when no unphysical cuts are 
considered, is given. If we denote by l)(s) a meson-meson partial wave with orbital angular 
momentum I the central result of ref.|| takes the form: 

Tl{s) = DiT) 
N(s) = 1 

i s ~ s i * Js thl - «)(«' - s o) 

Here p(s) = \q\/(8ir \fs), see eq. Q2,4| ), s t h r is the value of s for the threshold of the reaction, 
so is the subtraction point, a{so) a subtraction constant and each pole of the sum is referred 
to as a CDD pole after ref.0], with Sj the pole position and Ri the corresponding residue. 
Note that a CDD pole is a pole of the D(s) function so that it corresponds to a zero of the 
partial wave T)(s). 

In order to translate the result of eq.( |4.1| ) to the 7rN case, we first consider how the right- 
hand cut, given by the integral in eq.( |4.1| ), looks like in the complex VF-plane. Note that 
because W = \fs and since the square root is a double-valued function, we will have two 
cuts in the VF-plane. One of these cuts extends from W > W t hv = M + m to +oo and the 
other one from — oo to — Wthr- Therefore, the contour of integration in the complex s-plane 
used in ref.|| to obtain eq.([4.1|), Fig. 6, transforms in the W^-plane to the one shown in fig. 7. 
In this last figure, each line is denoted by two numbers i refers to the corresponding line 
in Fig. 6 and j = 1 means first sheet for W = y/s (ImW > 0) and j = 2 refers to the second 
one (ImVF <0). Then the dispersion integral in the W-plane for D(W), with a subtraction 
in Wo, will be: 

w ~ w ° { r dW > ^0^)11 r dw , Dm» 

2m \J Wthr (W'-W)(W'-Wo) J Wthi (W'-W)(W'-W ) 

+ r dW ' D{ ~ W ' hl r dW Di ~ W ' )l2 \ (4 2) 

Jwi* (W' + W)(W' + W ) J Wthi (W' + Wo)(W' + W)j { ' } 

where we have indicated by a subindex ij the line in Fig. 7 over which the integral has to 
be taken. Making use in the former expression of the fact that T(W) = T*(W*), where * 
means complex conjugation, we then have 



W- W Q 



7T 



dW 

w thT 



lmD(W') u a lmD(-W') 2 i 



{W> -W){W' -Wo) {W> + W){W + Wo) 



(4.3) 
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s-plane 



< j= 

(m+Mf 



Figure 6: Contour of integration for the Cauchy representation of D(s) eq. ( |4.1|) . 
The enclosing circle extending to infinity is not shown. 



Since the line denoted by 21 in Fig. 7 corresponds to the line 2 in fig. 6, then 

lmD(-W') 2 i = -lmD(W') n = p(W) 



(4.4) 



and hence we finally have 



W- W 



■>w iht 

"OO 



SO 



7T 



dW'lmD(W') u 
P(s') 



(W'-W){W'-W ) (W f + W)(W' + W ) 



so J 



(4.5) 



with so = Wq and s' = W 2 . Thus, we recover the same integral as in eq.(<Ll) but this time 
in the complex IF-plane. The remaining terms in eq. (|4.1|) are not altered: the subtraction 
constant is, of course, also present in the IF-variable case, and the CDD poles will be now 
in the form R^/iyV — Wj). Note that, as stated above, these poles represent just zeros of 
the Tf ± (W) partial wave and then have to appear as poles of the D(W) function but now 
in the VF-variable. Note that a complex CDD pole in the s-variable gives rise to four poles 
in the IF -variable. Hence, instead of eq. fl4.ip we will have for any 7rN partial wave T^{W): 



TL(W) 



R' 



^W-Wj 



+ a(s ) 



■so 



7T 



ds - 



p(s') 



*th, 



s)(s' - s ) 



(4.6) 



It is important to stress that in what follows we will fix the subtraction constant at a point 
so where it is real. This, however, is done for convenience and need not be done. For a(so) 
to be real, we have so < s t hri with s t h r the threshold of the reaction under consideration. 
Note also that eq. Q4.6D allows to accomodate the possible pole of Tf±(W) at IF = 0, due 
to the kinematical singularities Q, which give rise to a 1/IF factor in our definition of the 
partial wave amplitudes. 



A connection with the large N c limit [32] is very enlighting. In this limit of QCD the 
dominant terms of any meson-baryon amplitude are of order (D(N®) and given by the local 
terms and pole contributions. These are obtained in eq. fl4.6|) by the inverse of the CDD 
poles and the leading piece of the subtraction constant, which are of order O (N®). Note that 
the integral in this equation, which contains the unitarity loops in the s-channel, is of order 
0(N ( T 1 ) and hence is subleading in the large N c counting, as expected. As a consequence, 
we split the subtraction constant in two pieces: 



a(s ) =a L + a SL (s ) 



(4.7) 
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W-plane 
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(m+M) 



m+M 



Figure 7: Contour of integration for the Cauchy representation of D(W) eq. 
The enclosing circle extending to infinity is not shown. 



where a L is 0(N®) and independent of the subtraction point sq, while a («o) is 0(N ( 7 1 ) 
and is linked to the integral so that their sum is independent of the subtraction point so- 
Therefore, we introduce the function g(s) 



,SL 



(so) - - 

7T 



ds - 



p(s>) 



On the other hand, we also need the K function 



so J 



K(W) 



R' 



W- W! 

Hence, in the absence of unphysical cuts, we can write 
T l I ± =[K(W)- 1 + g(s)]- 1 



-i 



K{W) 



l + K(W)g(s) 



(4- 



(4.9) 



(4.10) 



in close analogy with a ET-matrix formalism. In this derivation we have considered the 
isospin limit. The extension of this is framework to include strong isospin violation is based 
on a coupled channel formalism as detailed in ref. |J. Further isospin breaking due to 
electromagnetic effects is harder to include and will not be considered in what follows. 



4.2 The left-hand cut 

We have discussed in detail the case of no unphysical cuts since we will take this case 
as our starting point. The unphysical cuts will be treated in a perturbative way. The 
perturbation will be done in a chiral loop expansion of the unphysical cuts. In this way, 
our final calculation will include their contributions up to one loop calculated at 0(p 3 ) in 
HBCHPTl^l 



This way of proceeding has already been used in ref.[|10| to study the strong 
WlWl scattering in the strongly interacting Higgs sector and is also being used in ref. [25] 
to study the Ktt, Krj S-wave scattering. 



Let us now discuss in detail our treatment of the left-hand cuts. 

any 7rN partial wave amplitude resulting from the contributions discussed in sect. 3. If we 



Denote by T^± +loop 



# 4 Strictly speaking the HBCHPT coun ting is referred to the g(W, t) and h(W, t) functions because 
the coefficient (E + m)/2 in front of eq. (|2.2| ) is not expanded. This is based on the fact that for a 
proper wavefunction renormalization and matching to the relativistic theory, one should keep this 
prefactor as it is given here. For a discussion on these issues, see ref. |33| . 
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assume the saturation of the HBCHPT counterterms up to 0(p 3 ) by resonance exchange 
and perform a non-relativistic expansion of ^ree+ioop' we wolU( l recover the HBCHPT 
amplitude up to 0(p 3 ) given in ref.Q- We now write 

tUW) = ^ 1 = N{W) (4-11) 
l±{ ' D(W) l + N(W)g(s) { ' 

with N(W) = K(W)(1 + 5N{W)). For the case of not including any unphysical contribu- 



tions, i.e. 5N(W) = 0, we recover eq.( 4.10 ). We expand this representation of the partial 
wave amplitudes in order to reproduce the one loop (D(p 3 ) HBCHPT result TJ. ree+lo : 

T^ +loop (W) = N(W) - N^Wfgis) , (4.12) 

where N\{W) is the 0(p) contribution of N(W) and is given by the lowest order HBCHPT 
amplitudes. Consequently 

N(W) = T^ +loop (W) + N^Wfgis) . (4.13) 

We note that T^ +loop satisfies perturbative unitarity up to 0(p 3 ) in the HBCHPT chiral 
counting, which can be traced back to the proper wave function renormalization of the tree 
graphs. However, to assure exact unitarity also for the loop graphs (which is of higher 
order in HBCHPT), we have to multiply g(s) in the previous equation by a correction 
factor, g= (2W)/(E + m) = 1 + 0{p 2 ). This assures a complete matching of the heavy 
baryon amplitude to the relativistic one with respect to the constraints from unitarity. 
Stated differently, if would not multiply the function g(s) by g, then the exact unitarity 



would be lost. Note that from the last equation and in eq.( |4.11| ), the N(W) and the D(W) 



functions satisfy the N/D equations up to the same order than our input T* Iee+li oop (W), 
that is, up to one loop calculated at 0(p 3 ) in HBCHPT. On the other hand, the D function 
satisfies them exactly for the right-hand cut. As a result, our final expression for any 
partial wave, eqs.(4.11) and ( [4.13 ), satisfies unitarity to all orders. We note that the low- 



energy representation could be further improved by matching to the fourth order HBCHPT 
amplitude, which should be available in the near future. 

TC(q) 

/ N 
/ \ 
/ \ 
/ \ 

L V 

N (P-q) 

Figure 8: Loop graph corresponding to the function g(s). 

Before finishing this section, we give an alternative expression for g(s) more suited than 
eq.( [4.8D for a direct comparison with HBCHPT. This new expression is derived making use 
of dimensional regularization. Since the function g(s) corresponds, up to a constant, to the 
loop graph with one pion and one nucleon propagator, see fig. 8, we can write: 

SL 1 f _ w, M 2 M 2 m 2 -M 2 + s, M 2 



g(s) = a^(A) + ? — ^ ^ -1 - - log + log -2- log — 

|<fl (m 2 + M 2 -s + 2W\q\\\ fA , A . 



W \m 2 +M 2 - s - 2W\q\ 
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In the former expression we have subtracted the piece 



w . M 2 

since this, together with the singular terms stemming from dimensional regularization: 

l < a » = t£{^4-> 4 '' +i " + 1 i} ' (4 ' 16> 

with d the number of space-time dimensions, is reabsorbed in the 0(p s ) HBCHPT coun- 
terterms |TJ. On the other hand, the subtraction constant a is fixed such that for a 
value of W below threshold and in its neighborhood, g(s) = 0. Note again that since 
we require d SL to be real, this matching has to be performed at or below the physical 
threshold, s t hr = (jn + M) 2 . It is known that the HBCHPT converges best in certain un- 
physical regions (but not in all), see e.g. ref.[[jl|. Therefore, in this energy domain we have 
Ty_ = T tIee+l = Ph^chpt + ^(p 4 )' which should be an accurate estimate of the true 
value of the partial wave since around threshold HBCHPT to 0(p 3 ) is convergent. We have 
fixed the regularization scale at A = m. Note that a change in the scale A is reabsorbed in 
a change of the constant d SL , such that at the end the value of g(s) is unaltered. 



5 Results 

In this section, we show and discuss the results. First, we explain how we have determined 
the parameters and which of these are strongly constrained by our fits. Then, we show some 
representative fits to the low-lying partial waves, exhibiting in particular the dependence 
on some of the coupling constants and mass parameters. 

5.1 Fitting procedure and parameters 

First, since we do not consider inelasticities, we restrict the fits to a cm. energy of W < 
1.3 GeV. We fit the two S- and the four P-waves, these are Su,Ssi and Pn, P13, P31 and 



P33, respectively. The fits are based on the Karlsruhe-Helsinki KA84 3£] and the VPI 
SM95 phase shift analyses. These PWAs do not give errors. One could assume an 
overall uncertainty of a few percent, but in that way the smaller partial waves would be too 
strongly weighted. Also, there should be some overall systematic error. To deal with these 
two effects and to give proper weight to the large partial waves (in particular the P33), we 
work with an error function for each partial wave of the form 



err (5) = ^ e 2 + 0.0025 5 2 , (5.1) 

with e s = or e s = y/2 and the phase shift 5 is given in radians. In this was a systematic 
error of e s and a statistical one of a 5% are added in quadrature. As will be shown in 
more detail below, the choice e s = gives more weight to the threshold region whereas for 
e s = y/2 we obtain a very precise description at higher energies, in particular for the A, and 
an overall better fit for the experimental data up to W = 1.3 GeV. Explicit figures for the 
fits are displayed in the next section. We now enumerate and discuss the parameters, which 
are essentially coupling constants and masses of the excited baryon and meson resonances. 
Fitting these parameters can be compared to a study in CHPT, where one fits the pertinent 
LECs. Clearly, for very low energies, this is analogous to assume resonance saturation of the 
LECs, which appears to be a sensible approach. However, due to the explicit propagators, a 
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fitting procedure as performed here, which is also unitary, can be extended to energies where 
HBCHPT can not work any more. We have performed various types of fits, leaving these 
parameters completely unlimited or prescribing limits on some of the parameters based 
on other phenomenological observations. These we will discuss in detail in what follows. 
Clearly, not limiting certain parameters can lead to rather large (unphysical) values for 
them. However, what we mostly can determine are products of couplings (with respect to 
some reference mass) and these normalized products always come out in very limited and 
sensible ranges. We should also stress that since our approach does not involve form factors, 
one has to be careful to compare our results to the ones obtained e.g. in meson-exchange 
models. Having made these general remarks, consider now the various fits. 

First, we discuss the generic results for the case e, = y/2. For the A(1232), we have the 
coupling constant and the parameter Z. We find that independently of how we perform 
the fits, Z ~ —0.1 and g n AN = 2.05. Note that the values for Z are well within the 



ranges of previous determinations, cf. refs.j34|, 35]. We note that the standard spin-isospin 
SU(4) using the Goldberger-Treiman relation, g n NN = OA'm/Fn, leads to (we normalize the 
coupling in terms of the charged pion mass) 

9nAN = ^=y- = 2 - 01 /M . (5.2) 

This means that the Delta parameters are well determined by the fit. We have also de- 
termined the complex pole in the energy plane corresponding to the P33 resonance and 
find 

(m A ,r A /2) = (1210,53) MeV , (5.3) 

in perfect agreement with dispersive analyses of pion-nucleon scattering [|| or pion pho- 
toproduction pOj . For the energy range we consider, the Roper only plays a role in the 
Pn partial wave. We have performed fits with and without Roper. Including the Roper, 
we find that for the coupling \fR, the fits prefer a somewhat larger value than the one 



used in ref.[21]. Next, we consider the scalar meson exchange. The pertinent parameters 
are Cd m ,Mg and c/g. For the constrained fits, we used c m = Cd = F n /^ ± 30% and 
Ms = 1200 MeV ±30%, with the central values discussed above. Although these parame- 
ters are not very well constrained by the fits, we observe a few remarkable features. First, 
the scalar mass prefers to come out large, which essentially means that the scalar-isoscalar 
exchange can be very well represented by a contact term. Also, the product of the scalar 
couplings to the the nucleon and the pions in units of the scalar mass squared turns out to 
be essentially constant, 

C -^- ~ (1.7... 2.1) • HT 3 MeV" 1 . (5.4) 



This number is about a factor of two larger than in ref.[ll] pointing towards a larger LEC 
c\ and consequently a larger pion-nucleon u-term. However, it is known that there are 
sizeable fourth order corrections to the scalar sector of (HB)CHPT and we therefore do not 
want to dwell on this point any further. At first sight, these findings are in disagreement 
with meson-exchange models or K-matrix approaches which like to have a light scalar- 
isoscalar meson. Note, however, that due to our chirally symmetric construction, the Smr 
coupling has a different momentum-dependence than what is usually assumed in these other 
approaches. To mimic this momentum-dependence, a light scalar is needed although the 
resulting fits are much worse in such case. For the p-meson, we have three parameters. 
The conventional ones are g p and K p . From these, only k p is reasonably well determined, 
K p = 6.1 ± 0.4 |}6|. There is a larger spread in the values obtained for g p either from the 
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width, the KFSR relation or fits to the nucleon electromagnetic form factors. To cope with 
this uncertainty and remembering again that we have no form factors at the vertices, in the 
constrained fits we limit g p between 5 and 8. For example, the standard Hohler-Pietarinen 



analysis 37] gives g p NN = 5.7. Note, however, that even in the fits without limits, g p 
always stays in this range. On the contrary, the value of k p is not well determined by the 
unconstrained fits. The additional parameter T+ is always negative and small. We introduce 
the dimensionless parameter T + = (M 2 Mp)/(v2 <7pAw) T + , and find —0.3 < T + < —0.1, 
i.e. the overall contribution of this additional term is small in most partial waves. In Sn and 
P13, however, varying T + in the range mentioned leads to visible changes in the phase shifts, 
in particular above W ~ 1.2 GeV. Therefore, it would be interesting to study the effect of 
this extra coupling in other approaches to 7rN scattering. The last parameter concerns the 
matching of the chiral amplitude to the unitarized amplitude at the subtraction point sq. 
Setting = m + aM, we always find that a = 0.16. 

We now turn to the fits with e s = 0. As already remarked, the A parameters are somewhat 
less well described, we get g-^AN = 2.25, Z between —0.2 and —0.1 and (mA,rA/2) = 



(1199,58) MeV for the complex pole. On the contrary, the coupling yR of the Roper is 



consistent with the one given in ref. [21|. The results for the scalar and the p meson 
parameters are similar to the ones for e s = \/2, again we note that the combination of 
the scalar couplings lies in the narrow range of (1.5 . . . 1.7) • 10 -3 MeV -1 . The matching 
of the chiral and unitarized amplitudes happens at larger values of a ~ 0.5. This is, of 
course, expected and leads to an improved description of the S-wave scattering lengths as 
compared to the case of e s = y/2. In the following section, we will tabulate some of the 
numbers discussed here and also display some representative fits. 

5.2 Representative fits 

After having discussed the genuine results in the previous paragraph, we give in table [j] 
the values of the various fitted quantities for four representative fits with e s = v2, based 
on the SM95 partial waves (using the KA84 analysis leads to very similar results). For all 
these fits, the x 2 /dof is very small. Although for the case of no limits on the parameters 
the x 2 /dof improves, these fits should only be considered illustrative since they lead to too 
large parameters in the scalar sector (we stress again that in such an approach as presented 
here one is only sensitive to products of couplings). The first fits are based on limiting the 
ranges allowed for the scalar couplings Cd, m , the scalar mass and the p couplings. In fits 3,4 
all these limits are lifted. Fits 2 and 4 differ from 1 and 3, respectively, in that the vector 
coupling Gy is reduced from 67 to 53 MeV (as explained before). The results are collected 
in table |l[ The x 2 /dof is of similar quality in all these fits. Note the remarkable fact that 
in all the fits, also for the case of the second table, the a parameter is always less than 1, 
as required by our final definition of the g(s) function, eq. ( |4.14 ). Note also that in the 



constrained fits, the product c m Cd comes out close to the large N c result of (38 MeV) 2 , cf. 
eq.( |3.14| ), although individually c m (cd) prefers somewhat larger (lower) values than F^/y/E. 
The scalar mass always tends to the upper end of the given range. For these four fits, we 
find c m g s /Ml = (2.1,1.9,1.7,1.7) • 10~ 3 MeV -1 , in order. If we, however, fix the scalar 
mass to 1200 MeV, Cd comes out closer to the large N c prediction and gs is smaller, but 
again we find that the aforementioned product of scalar couplings is 1.7 • 10 _3 MeV _1 . As 
noted before, our fits are only sensitive to the products Gyg p , Gyn p and GyT + . This is 
reflected by the fact that for fits 3 and 4, we find the same value for Gyg p and almost the 
same for fits 1 and 2. The corresponding S- and P-wave phase shifts are shown in fig.|9[ 
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Fit 1 


Fit 2 


Fit 3 


Fit 4 


9nAN 


2.05 


2.04 


2.05 


2.05 


z 


-0.16 


-0.08 


-0.05 


-0.05 


Vr 


0.67 


0.79 


0.79 


0.79 


f + 


-0.15 


-0.26 


-0.09 


-0.09 


c d [MeV] 


26.2 (25. ,50.) 


25.0 (25.,50.) 


5.9-10 5 


115.8 


c m [MeV] 


50.0 (25. ,50.) 


43.9 (25.,50.) 


1.0-10 6 


199.8 


M s [MeV] 


1560 (840,1560) 


1560 (840,1560) 


3.8-10 7 


63681. 


9s 


100.2+ 


104.3* 


2.4-10 6 + 


35159.+ 


G v [MeV] 


67* 


53* 


67 * 


53* 


9 P 


5.00 (5,8) 


5.00 (5,8) 


6.08 


7.70 


Kp 


5.70 (5.7,6.5) 


5.70 (5.7,6.5) 


3.53 


3.53 


a 


0.16 


0.16 


0.16 


0.16 


xVdof 


21.5/355 


15.6/355 


12.9/355 


12.0/355 



Table 1: Representative fits with e s = as explained in the text. The * denotes 
a fixed quantity and the + indicates that gs appears always multiplying c dm . For 
fits 1,2, the ranges on the various parameters as explained in the text are given in the 
round brackets. 





Fit 1 


Fit 2 


Fit 3 


Fit 4 




2.21 


2.25 


2.26 


2.26 


z 


-0.18 


-0.12 


-0.09 


-0.09 


Vr 


0.46 


0.56 


0.54 


0.54 


T + 


-0.19 


-0.23 


-0.08 


-0.08 


c d [MeV] 


25. (25. ,50.) 


25.7 (25. ,50.) 


5.7-10 5 


87847. 


c m [MeV] 


45.2 (25. ,50.) 


45.3 (25. ,50.) 


1.0-10 6 


126.7 


M s [MeV] 


1560 (840,1560) 


1560 (840,1560) 


3.9-10 7 


222.5 


9s 


92.6+ 


87.9+ 


2.3 -10 6t 


53682.+ 


G v [MeV] 


67* 


53* 


67 * 


53* 


9 P 


5.00 (5,8) 


5.62 (5,8) 


6.76 


8.55 




5.7 (5.7,6.5) 


5.7 (5.7,6.5) 


3.64 


3.64 


a 


0.65 


0.41 


0.51 


0.51 


X 2 /dof 


1.15 


1.00 


0.78 


0.78 



Table 2: Representative fits with e s = as explained in the text. The * denotes 
a fixed quantity and the + indicates that gs appears always multiplying c dm . For 
fits 1,2, the ranges on the various parameters as explained in the text are given in the 
round brackets. 
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Figure 9: Fits to phase shifts as a function of the cm. energy W for e s = \f2 and 
Gy = 53 MeV. The solid (dashed) lines refer to the constrained (unconstrained) fits 
as described in the text. 



We point out that in particular the P\\ and the P33 partial waves are well described, and 
a good overall description of data up to W = 1.3 GeV is achieved. If one calculates the 
S-wave scattering lengths and the P-wave scattering volumes, these come out in reasonable 
agreement with the data. 

We now turn to the fits with e s = 0. We have performed the same four fits as described 
above. The resulting parameters are collected in table |2[ These are similar to what was 
found before. The most striking differences are the Roper coupling, which comes out close to 
the value derived in ref.[^] and the larger value for the matching parameter a. Still, we find 
a < 1, as demanded by our definition of the matching function g(s). The resulting x 2 /dof 
is around 1 . Note that this number is only of limited meaning since it strongly depends on 
the choice of e s , which now is and for table [l] it was \/2, and on the assumed statistical 
uncertainty. The resulting S- and P-wave scattering lengths (volumes) are collected in 
table |3|. With the exception of aj~ + , they are in good agreement with the results based 
upon the Karlsruhe and VPI partial wave analysis. We stress that aj, is very sensitive 
to cancelations between contact term and loop contributions, which are individually much 



bigger than the total sum [|4TJ, and thus is only expected to be described precisely at fourth 
order. The resulting phase shifts for W < 1.3 GeV are shown in fig.|l(]. They are slightly 
worse than the ones for the fits with e s = y2. In particular, the P33 partial wave is not 
as precisely described as before, but still the complex pole is only slightly displaced as 
compared to the fits with e s = y/2. 
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Obs. 


Fit 1 


Fit 2 


Fit 3 


KA85 


SP98 


a d+ 


1.59 


1.13 


0.87 


-0.83 


0.0 ±0.1 


°o+ 


8.37 


8.31 


8.28 


9.17 


8.83 ±0.07 


af_ 


-4.92 


-4.85 


-5.03 


-5.53 


-5.33 ±0.17 


af + 


13.71 


13.87 


13.73 


13.27 


13.6 ±0.1 


a^[_ 


-1.23 


-1.27 


-1.21 


-1.13 


-1.00 ±0.10 


a i+ 


-8.08 


-8.15 


-8.10 


-8.13 


-7.47 ±0.13 



Table 3: Values of the S- and P-wave threshold parameters for the various fits as 
described in the text in comparison to the respective data (taken from ref. [0]). Fits 3 
and 4 give identical results within the precision displayed. Units are appropriate 
inverse powers of the pion mass times 10 -2 . 
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Figure 10: Fits to phase shifts as a function of the cm. energy W for e s = and 
GV — 53 MeV. The solid (dashed) lines refer to the constrained (unconstrained) fits 
as described in the text. 
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6 Conclusions 



In this paper, we have developed a novel chiral unitary approach for meson-baryon processes. 
As a first application, we have considered elastic pion-nucleon scattering from threshold up 
to the opening of the inelastic channels at W ~ 1.3 GeV. Making use of the time-honored 
N/D method, we have derived the most general structure of any pion-nucleon partial wave 
amplitude neglecting unphysical cuts. The latter are then treated in a perturbative chiral 
loop expansion based on the power counting of heavy baryon chiral perturbation theory. 
The most important and novel features and results of this approach can be summarized as 
follows: 

i) The approach is based on subtracted dispersion relations and is by construction rela- 
tivistic. No form factors, finite momentum cut-offs or regulator functions are needed 
to render the resummation finite. This is done by a subtraction of the dispersion 
relations below the physical threshold. 

ii) The amplitude is matched to the one of (heavy baryon) chiral perturbation theory 
at third order. This is done at an energy slightly below the physical threshold. In 
this region, the HBCHPT amplitude is expected to converge. It is straightforward to 
extend this matching to a fourth order and/or a relativistic CHPT amplitude. 

iii) The explicit resonance degrees of freedom contain free parameters (couplings and 
masses) which are determined from a best fit to the pion-nucleon partial waves. 
Some of these parameters, in particular the ones of the A, are very well determined. 
Couplings of the scalar and vector mesons can be pinned down less precisely. We have 
also pointed out a new coupling of the p not considered in conventional schemes. Our 
approach is consistent with the concept of resonance saturation of the low-energy 
constants. 

iv) The pion-nucleon phase shifts are well described below the onset of inelasticities. 
With the exception of the isoscalar S-wave scattering lengths, we also find a reason- 
able description of the 7rN threshold parameters. This can be improved by matching 
to a more precise (HB)CHPT amplitude. 

v) Our method is very general. It embodies (as a particular limiting case) the inverse 
amplitude method (IAM) |4^] recently applied to the 7rN scattering in |43|| . For a 
discussion about this point see refs. JlOj , [44[ . The chiral unitary approach pi], based 
on the lowest order tree level meson-baryon interactions, appears also as a special 



case 



vi) The extension of this approach to coupled channels (like irN — > r]N) and to the 
three-flavor sector (with coupled channels) is in principle straightforward. For the 
kaon-nucleon system, the matching to the chiral perturbation theory amplitude is 
more tricky due to the appearance of subthreshold resonances. Work along these 
lines is in progress. 
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